Laser machining simulation method, laser machining system having means for implementing the method, and computer program for implementing this method

ABSTRACT

A method for simulating laser machining of a material by a laser machining system comprising the following steps: providing a central unit with: information relating to the material to be machined: delta δ, threshold fluence, incubation coefficient S, complex refractive index n+ik, information about the laser processing system: information relating to a polarization, pulse energy E P , diameter of said machining laser beam at a focal point w, order of a Gaussian p, pulse repetition rate PRR n, wavelength; determining with said central unit on the basis of the information relating to said material to be machined and the laser machining system, a machining profile in two dimensions corresponding to the simulation of a machining of said material to be machined with said laser machining system.

TECHNICAL FIELD

According to a first aspect, the invention relates to a method for simulating laser machining of a material. According to a second aspect, the invention relates to a laser machining device for carrying out the machining of a material on the basis of a simulation obtained with the method of the first aspect of the invention.

BACKGROUND

Laser beams are used in particular for machining workpieces. It is possible to melt, evaporate or sublimate part of a material exposed to a laser beam. The machining of workpieces with laser beams enables highly efficient cutting, drilling, grooving, surface texturing for example.

Laser machining devices increasingly use laser beam pulses with a duration of less than one ns and a high energy per pulse of the order of J/cm². It is necessary to define and parameterize the laser machining device so as to obtain the desired target machining. It is then necessary to define the laser machining parameters for any newly defined target machining as well as for any new material.

When performing a new target machining and/or using a new material, it is necessary to perform preliminary machining experiments to set the laser machining parameters for achieving the target machining for the selected material. A large number of preliminary experiments is often necessary, which requires prolonged use of a laser machining device by an operator, as well as numerous material samples.

A technical problem that needs to be solved is to have a method for defining laser machining parameters without the use of a laser machining device as well as numerous material samples and their characterization. In particular, it would be useful to be able to predict the results of laser machining without having to use a laser machining system.

SUMMARY OF THE INVENTION

According to a first aspect, one of the purposes of the present invention is to provide a method for the reliable simulation of laser machining on the basis of laser machining parameters.

To this end, the inventors propose a method for simulating a laser machining of a material by a laser machining system capable of emitting a machining laser beam, said method comprising the following steps:

a) providing a central unit with:

-   -   information relating to the material to be machined comprising:     -   a delta δ relative to a penetration depth of the machining laser         beam in the material,     -   a threshold fluence F_(th) relative to a minimum energy density         of the laser beam allowing machining,     -   an incubation coefficient S, between 0 and 1, relating to a         dependence of the threshold fluence F_(th) on a number of pulses         of the laser beam at a single point,     -   a complex refractive index n+ik,     -   information relating to said laser machining system comprising:     -   information relating to a polarization of said machining laser         beam,     -   a pulse energy of said machining laser beam E_(P),     -   a diameter of said machining laser beam at a focal point w,     -   an order of a Gaussian p of said machining laser beam between 1         and 20,     -   a pulse repetition rate PRR of pulses n of the machining laser         beam,     -   a wavelength of the machining laser beam;

b) determining with said central unit on the basis of:

-   -   said information relating to said material to be machined, and,     -   said information relating to said laser machining system,

a machining profile in two and/or three dimensions corresponding to the simulation of a machining of said material to be machined with said laser machining system.

The simulation method according to the invention makes it possible to solve the technical problem thanks, in particular, to the proper consideration of laser-material interaction, and to its implementation in said simulation method. The simulation method of the invention proposes a numerical approach for the simulation of complex pysitic phenomena so as to obtain simulation results very close to the experimental results obtained with the same laser machining parameters. Indeed, the simulation tools of the prior art use an analytical approach so that the resolution of the analytical models is very often very demanding in terms of computing capacity and does not allow certain simulations, for example when a simulation presents an incident laser beam with an incidence that is not normal to the surface of the material.

The method of the invention allows the determination of optimal machining parameters and preferably of optimal machining parameters allowing a piece to be machined according to a target machining as close as possible. Such a method makes it possible to guarantee good machining quality without having to resort to numerous and costly machining parameter adjustments. Indeed, the method proposed by the invention makes it possible to take into account the physical parameters of interaction between the laser beam and the material to be machined during machining, during previous machining experiments or from a database containing information relating to a laser source or a material to be machined. Preferably, the method of the invention foresees to provide before step a), a step of providing a central unit.

The inventors have observed in a surprising manner that the taking into account by the control unit of the four material parameters in relation to laser machining (delta, threshold fluence, incubation coefficient, complex refractive index) makes it possible to obtain machining simulations that are very close to the machining to be achieved without having recourse to empirical experiments and for extended ranges of use of a laser machining device.

Preferentially, the information relating to the material to be machined is information relating to the interaction of the laser machining beam with said material to be machined. Preferentially, the information relating to the polarization of said machining laser beam selected from a choice of vertical or horizontal P/S polarization polarization. Preferentially, the pulse energy of the machining laser beam Ep is expressed in μJ. Preferentially, the diameter of the machining laser beam at a focal point w corresponds to the waist size w. Preferentially, the order of one Gaussian of said machining laser beam is between 1 and 10.

According to a preferred embodiment of the invention, the material parameters are stored in the database after an optimization step. The step of optimizing the material parameters in the database includes a comparison of experimentally measured material parameters and material parameters estimated on the basis of an analytical model. The comparison based on the analysis of the R², which must be between 0.92 and 1, allows an optimization of a material parameter in the database so as to minimize errors during simulation due to the reduction of possible errors on the material parameters recorded in the database. In particular, the step of determining the material parameter delta δ is done by performing a comparison between the experimental measurements and the model of the ablation depth as a function of fluence. The value of the delta δ is determined after several iterations until the superposition of these two curves measured by the R² is between 0.92 and 1.

The numerical approach of the invention makes it possible to successively simulate the arrival of each pulse of the machining laser beam on the surface of the material. The topology of the material surface being updated before the next pulse. Thus the simulation of the invention is based on an understanding of the laser-material interaction of each pulse, which makes it particularly accurate and realistic in comparison with other prior art analytical methods.

The method of the invention enables the laser machining parameters to be predicted by simulation and allows them to be applied to the workpiece by means of a laser machining device, in order to obtain a target type of machining (engraving, cutting, drilling, welding). These characteristics/properties of each material (both physical and learned functions) make it possible to predict the optimal machining parameters (laser conditions) to be applied to obtain a desired machining.

Preferably,

the laser machining system is a system capable of emitting a laser machining beam moving in a direction x substantially parallel to a machining surface of said material;

step a) further comprises providing the following information relating to said laser machining system:

-   -   a distance of the focus point of the laser machining beam from         the surface of the material to be machined,     -   a speed (v) of movement of the laser machining beam relative to         said material to be machined     -   an angle of incidence of the machining laser beam with respect         to said surface of the material to be machined     -   a number of lines to be machined     -   a distance between said lines previously defined     -   number of passes of the machining laser beam on each line to be         machined.

For example, the distance of the focus point of the laser machining beam from the surface of the material to be machined is 0 when the focus point is on the surface. Preferably, the moving speed is set in the direction x. For example, the angle of incidence of the machining laser beam with respect to the surface of the material to be machined is 0° when the laser beam makes an angle of 90° with the surface of the substrate. For example, a number of lines to be machined is at least 1 and can tend towards infinity. Preferably, the number of lines to be machined allows a trench width to be defined in relation to the distance parameter between the lines. Preferentially, the distance between said lines or pitch is a distance from the center of one line to the center of another adjacent line. For example, the number of passes of the machining laser beam on each line to be machined is a number of layers. A number of passes of the laser beam equal to 10 means that the laser beam passes 10 times along the same line.

Preferably, the laser machining system is a system capable of emitting a laser machining beam in rotation about an axis of rotation substantially perpendicular to a machining surface of said material and in that, the step a) further comprises providing information on said laser machining system:

-   -   a rotational speed ω of said laser beam about said axis of         rotation,     -   a distance BFG from a surface of said material to be machined         for which the rotating laser beam describes a fixed spot for all         rotating laser machining beam positions,     -   a distance BFI from a surface of said material to be machined at         which the rotating laser beam is focused,     -   an angle of incidence β of the rotating machining laser beam for         all positions of the rotating laser beam with respect to the         normal to the surface of the material     -   an activation time of the rotating machining laser beam.

Preferably, the rotational speed of the laser machining beam is the rotational speed of an optical element allowing the rotation of the laser machining beam. For example, such an optical element is a mirror, a prism, etc. The distance BFG from a surface of said material to be machined for which the rotating laser beam describes a fixed spot for all positions of the rotating laser machining beam, i.e. a point around which the laser machining beam rotates. Preferably a distance BFI from a surface of said material to be machined at which the rotating laser beam is focused, e.g. a distance BFI equal to 0 means that the laser beam is focused on the surface of the material for all positions of the rotating laser beam. Preferably, the distance BFI from a flat surface of a material to be machined is a plane. Preferably, an angle of incidence β of the rotating machining laser beam for all positions of the rotating laser beam relative to the normal to surface of the material, 0° relative to the surface means that the rotating laser beam describes an angle of 90° from the surface or 0° from the normal to the surface of the material. Preferably, an activation time of the rotating machining laser beam depends on the thickness and type of material to be machined and other laser parameters. This embodiment is particularly advantageous because it allows to simulate with precision the laser parameters which will enable to have straight or negative flanks (the orientation of the flanks in general).

Preferably, the machining laser beam is a machining laser beam in rotation about an axis of rotation and a point of rotation located at said distance BFG from said surface of said material to be machined.

Preferably, that step b) of determining said machining profile in two and/or three dimensions comprises a step of determining a precession radius r_(P) of a machining of said material by said rotating machining laser beam by the following formula:

$r_{P} = {\left( {\frac{BFG}{\cos \beta} - z} \right)\tan {\beta.}}$

Preferably two successive pulses n of the rotating machining beam are separated by a distance dx, preferably along a circle of radius r_(P).

Preferably, step b) of determining said machining profile in two dimensions comprises a step of determining an ablated crater radius r_(c) for a pulse of the machining laser beam by the following formula:

$r_{c} = {w_{0}{\sqrt{\ln \left( \frac{F_{0}}{F_{th}} \right)}.}}$

Preferably, the step b) of determining said machining profile in two dimensions comprises a step of calculating an ablation depth z_(n) by the following formula:

$z_{n} = {- {\delta \left\lbrack {{\ln \left( \frac{F_{0}}{F_{th}\left( {n + {N{R/2}}} \right)} \right)} - {2\frac{n^{2}dx^{2}}{w^{2}}} - {2\left( \frac{\left( {y + y_{c}} \right)^{2}}{w^{2}} \right)^{p}}} \right\rbrack}}$

Preferably, this machining profile is defined in two dimensions according to:

an axis y representing a direction substantially parallel to the surface of the material to be machined, and,

an axis z representing a direction substantially perpendicular to the surface of the material to be machined, the axis z corresponding to an ablation depth z_(n) with respect to the surface of the material.

Preferably, the delta δ is a constant parameter. The use of a constant parameter for delta δ is particularly advantageous because it allows to have a constant delta parameter δ for each material, regardless of the machining conditions applied. The inventors found that the use of a constant delta δ was contrary to the teachings of a large number of scientific publications. In particular, it has often been accepted that the delta δ of a material should vary as a function of the number of pulses of a laser beam, whereas the delta δ is an intrinsic property of the material and should not be considered as a fitting variable of a model.

Preferably, the central unit is configured to execute a computer program for the determination of said machining profile in two and/or three dimensions.

Preferably, the database of the invention allows the storage of at least one type of information as defined below: material characteristics, response of a material to one or more laser pulses. For example, a database includes material characteristics.

Preferably, said information relating to said material to be machined is taken from a material database, wherein each material comprises the following information:

-   -   a delta δ relative to a penetration depth of the machining laser         beam into the material,     -   a threshold fluence F_(th) relative to a minimum energy density         of the laser beam allowing a machining,     -   an incubation coefficient S between 0 and 1, relating to the         dependence of the threshold fluence on the number of pulses of         the laser beam at the same point,     -   a complex refractive index n+ik,

and in that said central unit comprises communication means for communicating with said material database.

Preferably the method of the invention enables and or applies to carry out simulations for observing a machining obtainable with parameters defined by a user. This is particularly used in order to be able to ascertain the feasibility of a machining or to define machining sequences.

The central unit is preferably set up to run a computer program for simulating a laser machining profile in two dimensions.

Preferably the method of the invention includes a step of providing the central unit with a computer program for executing modeling means and/or simulation means. Preferably, the modelling/simulation and/or machine learning means comprise a model based on algorithms for relating physical parameters to each other. An example of a computer program is software. Preferably, the modeling and/or machine learning means are used to generate data that can be stored in the database.

Preferably, the central unit configured to execute said computer program is accessible on dematerialized storage and computing means, e.g. a cloud, in particular a platform as a service and more particularly a software as a service.

The software as a service provides that this computer program is made available to the user. The computer program can be manipulated using a web browser or installed on a leasehold way on a central unit. The user does not have to worry about making updates, adding security patches and ensuring availability of the service.

For example, the method according to the invention is an internet platform to help a user who wants to perform machining with a laser machining device to choose the right machining parameters. Preferably the method according to the invention is a decision support tool. Preferably, the machining parameters proposed by the user to the central unit during the simulation are directly transmitted to the control unit after validation by the user.

Preferably, said central unit comprises modeling means and simulation means so as to be able to determine said machining profile in two dimensions on the basis of information relating to said material to be machined from said database, and on the basis of information relating to said laser machining system.

Preferably, the method further comprises the following additional steps:

a) providing a laser machining system comprising:

-   -   said material to be machined;     -   said laser machining system comprising:     -   a laser source for emitting a laser machining beam onto said         material to be machined;     -   a control unit for controlling the emission of said laser         machining beam from said laser source;     -   said central unit for controlling said control unit;     -   said communication means for communicating with the database;     -   communication means for communication from said central unit to         said database;

b) machining said material with said laser source configured with said laser machining parameters.

Preferably the laser machining device and the analysis unit in particular comprises sensors installed at many points of the machine for measuring the power in the machine and at the target, a beam analyser, a wavefront analyser (e.g. with a Shack-Hartmann), an M² measuring instrument, an instrument for calculating beam propagation with prediction of the position and size of the waists (position where the beam width is minimal) along the beam propagation, means for corroborating a measurement of the size of the spot and position on the workpiece (e.g. with a Nanoscan). Preferably the laser source is an ultra-short laser source, i.e. for sending laser pulses of less than 100 ps.

Preferably, the machining system includes means for analysing the result in real time, e.g. with an OCT (Optical Coherent Tomography), or after machining with an interferometer. The means for analysing the result in real time, for example with OCT and/or confocal OCT, whose different measured dimensions allow us to deduce the physical characteristics of the materials, and also to check the result obtained and to deduce an error, according to which we can apply a correction of one or more parameters calculated by the algorithm used for the simulation.

The method of the invention presented in this particular embodiment also aims at providing and feeding a material database which then allows the prediction by simulation of the optimal laser parameters specific to each material, according to criteria or within error intervals of the desired result. The invention includes a laser machining system which can be equipped with different laser sources (laser system apparatus), a laser beam management system, a machining unit allowing the relative displacement of the laser beam with respect to the workpiece and an analysis unit allowing a result to be obtained with respect to a machining carried out. An advantage of the invention is to have a material database allowing the storage of material characteristics. This material database makes it possible to have material characteristics from the start of a new machining in order to benefit from the results of past experiments when these are compatible with the current experiment, in particular when the material to be machined is known in the database.

Preferably, said laser machining device further comprises:

an analysis unit of the state of said material to be machined connected to said central unit to provide a machining result,

and that the following additional steps are implemented:

c) acquiring with said analysis unit after step d), machining results;

d) transmitting said machining results and said machining parameters to said central unit, said central unit being configured to determine information relating to the material to be machined;

e) enriching said material database with said information relating to the material to be machined from said central unit.

This preferred embodiment of the invention, allows an enrichment of the material database according to the invention for which an optimization algorithm/an optimization as described can be applied by the central unit so as to reduce errors due to measurement or other external factors.

Data Used by the Central Unit

The data used by the central processing unit shall comprise at least one data type with the following characteristic: data describing a target result or a result to be achieved, user parameters, machine-related characteristics, instructions corresponding to predefined tests, parameters related to a material to be machined, parameters related to a material stored in a database, parameters related to a material stored in a database and obtained during a previous machining, algorithms comprising learning functions, optimal machining parameters, data from means for measuring the laser beam, data from means for analysing the result. Preferably the types of data described above are annotated.

Data Cleaning

A cleaning step of the data collected or present in the database is provided for according to a preferred embodiment of the invention. This step ensures that data collected from a database or measured online are consistent, with no outliers or missing values.

Target and Predefined Results

The specific task to be performed or result to be achieved or target machining result or predefined target result according to the invention corresponds to the problem that a user of the present invention seeks to solve by modelling/simulating the phenomenon. The invention makes it possible to consider a plurality of target machining results, each of which may require the choice of different simulation algorithms.

Preferably the central unit of the laser machining system includes the control of all the units included in it and in particular at least: the laser source, the beam management device, the beam movement/focusing unit, the analysis unit, a database, a material database, an observation unit, a unit, a unit for implementing modeling/simulation, a unit including predefined test characteristics. A central unit is for example a computer with several input and output interfaces. The central unit is preferably connected to a network, for example to the Internet, allowing for example the relocation of units to implement modeling/simulation. The advantage that the units for carrying out the modeling/simulation are not physically present in the laser machining system is that a central unit is available which does not require the hardware resources, e.g. computing power required to carry out the modeling/simulation. This also allows to have a laser machining system with a central unit dedicated to laser machining. Preferably the database(s) are also stored on a central server, preferably the central server being located in the same place as the units to implement the modeling/simulation.

Preferably said database includes information related to said material of at least one of the following types: learned function, physical characteristics of a material.

Preferably, said physical characteristics of a material include at least one of the following characteristics: a delta, a threshold fluence, a maximum efficiency fluence.

Preferably, said analysis unit of the state of said material to be machined is an optical detection unit.

Preferably the laser machining parameters include at least one of the following parameters:

a scanning speed of the laser beam over the material to be machined;

a number of pulses per second emitted by said laser source;

a wavelength emitted by the laser source;

a spot size of the laser source on the material to be machined;

an energy per pulse emitted by said laser source;

angle of incidence of the beam on the workpiece;

a number of passes or layers to be machined;

the distance between two adjacent joint lines;

number of revolutions per minute of the optical rotation system.

Preferably, the control unit is further configured to control said light beam from said laser source to said material to be machined.

Preferably, the method further comprises the following additional step:

f) acquiring information relating to said material to be machined enabling identification of said material to be machined by the central unit.

Preferably, the method further comprises the following additional step:

g) further determining a machining strategy on the basis of said machining of said material to be machined.

Preferably, the method further comprises the following additional step:

h) irradiating said material to be machined with said light beam with predetermined light beam parameters;

i) acquiring by the optical detection unit a result generated by the irradiation of said material to be machined with said light beam;

j) transmitting said result to the central unit and accessing said result generated by irradiating said material to be machined with said light beam;

k) extracting from said result by means of the central unit at least one laser machining parameter relating to a material;

l) recording said at least one laser machining parameter relating to a material in said database.

The inventors propose, according to a second aspect, a laser machining system comprising means adapted to perform the steps of the method according to the first aspect of the invention.

The different variants and advantages described for the method according to the first aspect of the invention apply to the system of the second aspect, mutadis mutandis.

Preferably, the system further comprises an optical unit for directing said laser beam towards said material to be machined.

Preferably, said optical unit for directing said laser beam towards said material to be machined allows precession of said laser beam.

The inventors propose, according to a third aspect, a computer program for implementing the method according to the first aspect of the invention.

The different variants and advantages described for the method according to the first aspect of the invention and for the system according to the second aspect of the invention apply to the computer program of the third aspect, mutadis mutandis.

Preferably, the computer program comprising the instructions which cause the laser machining system according to the second aspect of the invention to perform the steps of the method according to the first aspect of the invention.

The inventors propose, according to a fourth aspect, a computer-readable medium on which the computer program according to the third aspect of the invention is recorded.

The different variants and advantages described for the method according to the first aspect of the invention, for the system according to the second aspect of the invention and for the computer program according to the third aspect of the invention apply to the computer-readable medium of the fourth aspect, mutadis mutandis.

Determining the Laser Machining Parameters on the Basis of the Laser Machining Parameters in Relation to the Material of a Workpiece

When a decision is made to machine a workpiece, a target machining of this workpiece is set. Preferably a target machining of a workpiece is a hole, a groove, a channel, a cutout or other machining obtained by ablation, through or non-through machining that the operator wishes to perform in a workpiece. Preferably, a target machining is defined in a two-dimensional space, and more preferably in a three-dimensional space. For example, a target machining defines a volume of material to be ablated. A target machining may, for example, correspond to a cut or a change in the optical properties of a transparent material.

In order to perform the target machining with the smallest possible machining tolerances, the machining parameters must be determined according to the nature of the workpiece, i.e. its composition, original shape and the properties of the laser source used. The method of the invention allows the determination of laser machining parameters without necessarily resorting to empirical tests beforehand. In fact, thanks to the use of the different models which will be detailed later and by knowing the interaction parameters of the material of the workpiece when it is irradiated by the light beam, it is possible to determine the laser machining parameters. The laser machining parameters determined when implemented according to the method of the invention make it possible to obtain the machining which would be obtained when machining with a laser machining system.

The laser machining parameters are preferably defined for a laser machining device and include but are not limited to:

a scanning speed of the laser beam over the workpiece;

a number of pulses per second emitted by said laser source;

a wavelength emitted by the laser source;

a spot size of the laser source on the material to be machined;

an amount of energy per pulse emitted by said laser source;

a pulse duration emitted by said laser source.

A workpiece is defined by the material of which it is made and its dimensions. For example, the dimensions can be defined in a multidimensional space in order to clearly define the material locations and the locations where no material is present. A dimension can for example define a roughness of the workpiece, for example a roughness can be defined for each of the faces of the workpiece when it is defined by faces.

It is possible to determine the δ from the equations described below, however it is not always easy to have all the parameters, constants or coefficients necessary to make this calculation. Thanks to the method of the invention it is possible to determine a 6 experimentally while taking into account the variability of certain parameters. To do this, scans must be produced over the surface of the material, with constant power values and increasing speed values, and with constant speed values and increasing power values, as shown in FIGS. 7a and 7b by way of example. Preferably, the maximum depths should be measured and compared with the values obtained with the corresponding model described below. The value δ used in the model must be adjusted until a good coincidence between the experimental points and the modeled points is obtained.

More generally, the coefficient δ, described as the absorption coefficient and used mathematically as a calibration factor. It is known from the prior art that the coefficient δ varies with the number of accumulated pulses due to an incubation coefficient S, so in the same way as for the ablation threshold, the coefficient δ is: δ(N)=δ(1)N^((S−2)).

In the numerical simulation of the invention, the coefficient δ is considered constant throughout the numerical calculation of the ablation profile, which gives results quite similar to those obtained experimentally with the same processing parameters.

Determination of the Ablation Threshold Using the D² Method.

Ablation threshold determination by the D² method is the most accepted and widely used method for determining the ablation threshold of all types of materials. This method is used for Gaussian beam energy profiles (see FIG. 1). This method has the advantage of offering the possibility to estimate the radius of the Gaussian beam co at the surface of the material and the incubation coefficient.

The fluence (or energy density) in a plane perpendicular to the laser beam axis is given by:

$\begin{matrix} {{{F(r)} = {F_{0}{\exp\left( \frac{{- 2}r^{2}}{\omega^{2}} \right)}}},} & (1) \end{matrix}$

where r is the distance to the beam axis, ω² is the radius of the beam at the surface of the material at e⁻² of the maximum intensity. If the beam is focused on the surface, ω=ω₀. F₀ is the maximum fluence (at r=0) and is calculated and given by:

$\begin{matrix} {F_{0} = {\frac{2E_{P}}{\pi \omega^{2}}..}} & (2) \end{matrix}$

The diameter of an ablation crater produced with one or more laser pulses can be described as a function of the maximum fluence:

$\begin{matrix} {{D^{2} = {2\omega^{2}{\ln \left( \frac{F_{0}}{F_{th}} \right)}}},} & (3) \end{matrix}$

where F_(th) is the ablation threshold. Since F₀ increases linearly with E_(P) the radius of the co beam can be estimated from the graph of crater diameters D² as a function of the logarithm of the pulse energy, using equation 3. The radius co can be calculated from the slope of a curve defined by the experimental results:

$\begin{matrix} {{\omega = \sqrt{\frac{slope}{2}}}.} & (4) \end{matrix}$

The ablation threshold can be calculated by extrapolating the craters diameter at power 2 to a crater diameter equal to zero in a graph representing the crater diameter at power 2 (D²) as a function of the logarithm of the pulse energy.

In general, the ablation threshold is reduced by increasing the number of pulses incident on the same point on the surface of the material to be machined. This behaviour is attributed to the increase in radiation absorptivity caused by the accumulation of defects created with each pulse, even at fluences lower than the threshold fluence. Such an increase in radiation absorptivity can be characterized by an incubation factor or coefficient, expressed by the following equations (5) and (6):

F _(th)(N)=F _(th)(1)×N ^(S−1)  (5)

When F_(th)(N) is the fluence of the ablation threshold for N pulses, F_(th)(1) is the fluence of the ablation threshold for a single pulse and S is the incubation coefficient. S=1 means that there is no incubation, i.e. the ablation threshold is independent of the number of pulses. For most materials S<1, the ablation threshold decreases with increasing number of pulses.

In general, the incubation or the decrease in threshold fluence is more pronounced during the first one hundred pulses for a large majority of materials and particularly for the most commonly used materials.

To apply the D² method to find the ablation threshold, the incubation coefficient and, if necessary, the beam radius, a crater matrix with increasing pulse numbers and pulse energies must be produced as shown in FIG. 2. In order to have sufficient experimental points, at least 20 energy values, and several pulse number values less than 100 pulses are required. The diameter of each crater must be measured from images obtained with an optical or electron microscope.

Determination of the Ablation Threshold Using the Diagonal Scan Method.

The Diagonal-Scan or D-Scan method is a geometric method for determining the value of the ablation threshold. The D-Scan method is an alternative method to the D² method and has the advantage of being more efficient and faster in terms of experiment time and number of measurements to be performed. The disadvantage of this method compared to D² method is related to the method for determining the beam radius, which is always important for measurements and calculation of the parameters.

The Diagonal-Scan method involves moving the sample diagonally in relation to the focal point of the beam using different beam energies and a variable number of pulses, while changing the speed of movement (FIG. 3). The experiments enabling to implement the D-Scan method involve the determination of the ablated surface. The ablated surface corresponds to the points where the radiation power has exceeded a critical value defined by:

P _(crit)=½eπω ₀ ² I _(th),  (7)

where I_(th) is the ablation threshold in terms of intensity, given by:

$\begin{matrix} {{I_{th} = \frac{P_{0}}{e\pi \rho_{{ma}\; x}^{2}}},} & (8) \end{matrix}$

where P₀ is the maximum pulse power and ρ_(max) is the maximum width of the ablated region (FIG. 3). The ablation threshold is given by:

F _(th) =I _(th)τ_(p),  (10)

Where τ_(p) is the pulse duration.

To apply the D-Scan method to find the ablation threshold and the incubation coefficient, scans must be produced as shown in FIG. 3, with different energy values and speeds of movement in order to have measurements for different numbers of superimposed pulses. The number of pulses superimposed on the points corresponding to the wider part of the lobes is given by:

$\begin{matrix} {{N = {\sqrt{\pi}\frac{f \times \rho_{{ma}\; x}}{v_{y}}}},} & (11) \end{matrix}$

where f is the pulse frequency and v_(y) the speed in the direction of the axis y.

Determination of δ Using the Skin Depth Method

During the interaction between an ultra-short laser pulse and a solid, the electrons do not have time to transfer their energies to the ions, because the electron-phonon coupling time is longer than the pulse duration. Under these conditions, the electron density can be assumed to remain constant and the electromagnetic field in the solid can be calculated using Maxwell's equations in combination with the material equations. It can also be considered that in this case, the electric field of the laser on the material surface decreases exponentially with the depth in relation to the material surface:

$\begin{matrix} {{{E(x)} = {{E(0)}{\exp \left( {- \frac{x}{\delta}} \right)}}},} & (12) \end{matrix}$

where —6 is the penetration depth or the skin depth. The surface of the material corresponds to x=0 and the equation (12) is valid for x>0. In general, the absorption depth is expressed as:

$\begin{matrix} {{\delta = \frac{c}{\omega_{L}\kappa}},} & (13) \end{matrix}$

Where κ is the imaginary part of the refractive index,

N=√{square root over (ε)}=n+iκ,  (15)

where ε=ε′+iε″ is the dielectric function and ω_(L) is the frequency of the laser. The dielectric function can be considered in the Drude formula for the following calculations:

$\begin{matrix} {{ɛ = {1 - \frac{\omega_{pe}^{2}}{{\omega \left( {\omega + {iv}_{eff}} \right)}^{\prime}}}}{where}} & (16) \\ {\omega_{pe} = \sqrt{\frac{4\; \pi \; e^{2}n_{e}}{m_{e}}}} & (17) \end{matrix}$

ω_(pe) is the frequency of the electron plasma, v_(eff) is the effective collision frequency between the electrons and the lattice of ions, and m_(e) and n_(e) are the mass and density of the conduction electrons, respectively. In the case of a high collision rate, v_(eff)>>ω and thus ε″>>ε′, the equation (13) can be reduced to the following form:

$\begin{matrix} {{\delta = {\frac{c}{\omega_{L}\kappa} \approx {\frac{c}{\omega_{pe}}\sqrt{\frac{2v_{eff}}{\omega_{L}}}}}}.} & (18) \end{matrix}$

Determination of δ by Ablation Depth

Under non-equilibrium conditions for the determination of δ by the skin depth method, the electron thermal conduction time (t_(heat)), or the time for the electrons to reach thermal equilibrium in the region in the skin layer is also longer than the laser pulse duration. The time t_(heat) can be approximated by the following equation, with an expression of the thermal diffusion coefficient as follows:

$\begin{matrix} {{{t_{heat} \approx \frac{\delta^{2}}{\kappa_{T}}};}{{\kappa = \frac{l_{e}v_{e}}{3}},}} & (19) \end{matrix}$

where K_(T) is the thermal diffusion coefficient, and l_(e) and v_(e) are the mean free electron path and its speed, respectively.

The laws of energy conservation in this case have a form that describes the variation in electron temperature caused by the absorption of energy in a skin layer:

$\begin{matrix} {{{{c_{e}\left( T_{e} \right)}n_{e}\frac{\partial T_{e}}{\partial t}} = {- \frac{\partial Q}{\partial x}}},} & (20) \end{matrix}$

where Q is the energy flow absorbed in the skin layer:

$\begin{matrix} {{Q = {Al_{0}{\exp \left( {- \frac{2x}{\delta}} \right)}}},} & (21) \end{matrix}$

where A=I/I₀ is the absorption coefficient, I₀=cE²/4π is the incident laser intensity. The absorption coefficient and the skin depth are therefore a function of the laser frequency ω_(L) of the conduction electron density n_(e), the effective frequency of electron-electron and electron-ion collision v_(eff), of the angle of incidence and of the polarization of the laser.

The depth of a crater z produced with an ultrafast laser with fluence F is of the order of the penetration depth of the beam, according to the equations:

$\begin{matrix} {{z = {\delta \; \ln \; \frac{F}{F_{th}}}},} & (22) \end{matrix}$

where F_(th) is a threshold fluence for the ablation and δ is a length having several interpretations depending on the approximations in the two temperature equations and the fluence range. Typically, in the low fluence domain, δ is related to the optical absorption coefficient and, for higher fluences, to thermal parameters such as electron scattering and the time interaction duration. For many materials, there are different values of F_(th) and δ at high and low fluences, with the boundary depending on the material. Each material is then represented by a set of parameters (δ, F_(th)) in a given fluence zone. The laser fluence F for a laser beam with a Gaussian intensity profile is given by:

$\begin{matrix} {{F(r)} = {F_{0}\exp \; \left( {- \frac{2r^{2}}{\omega_{0}^{2}}} \right)}} & (23) \end{matrix}$

where r is the radial distance (distance from the center of the laser beam), F₀ is the fluence peak or the fluence at the center of the beam, and ω₀ is the radius of the beam at 1/e² of the maximum intensity. The maximum fluence F₀ can be obtained from the pulse energy E_(p):

$\begin{matrix} {{F_{0} = \frac{2E_{p}}{\pi \omega_{0}^{2}}}.} & (24) \end{matrix}$

The depth of pulse ablation after irradiation with a Gaussian beam is given by the relation:

$\begin{matrix} {{{z_{G}(r)} = {{\delta {\ln \left( \frac{F_{0}}{F_{th}} \right)}} - {2\frac{r^{2}}{\omega_{0}^{2}}}}}.} & (25) \end{matrix}$

This model is based solely on the response of the material to the applied fluence. All geometric effects described by this approach will be related to the temporal or spatial distribution of fluence.

In the case of multi-pulse irradiation, it must be taken into account that the surface topography is no longer flat but becomes paraboloid, its depth increasing after each consecutive laser pulse. Therefore, a few considerations are essential to achieve a better prediction of the ablation dimensions. The first condition that should not be overlooked is the incubation, or lowering the ablation threshold with the number of pulses. Secondly, the effects of local surface angle on energy absorption and reflectivity must be taken into account.

Finally, and since multi-pulse irradiation is typically used to obtain deeper craters or laser scans, the variation of the Gaussian beam of the fluence with the depth must also be considered.

This description is valid for relatively low fluences with respect to the ablation threshold. For higher fluences, thermal parameters such as electron scattering or the duration of the laser-matter interaction must also be considered.

Line Ablation

In most applications, a static description is not relevant if the sample or beam is moving. A numerical description can be used for dynamic situations. Consider a motion in the direction x with a speed v. Since the pulse repetition rate of the laser is PRR, the distance Δx between two successive pulses is:

$\begin{matrix} {{{\Delta x} = \frac{v}{PRR}}.} & (26) \end{matrix}$

In this case, it is not possible to assume a radial symmetry and the fluence is then a function of the position in x and in y. It is then necessary to enter a number of pulses N as seen by the material for one pass. The ablation depth can then be expressed according to the following equation:

$\begin{matrix} {\mspace{85mu} {{{z_{k}\left( {0,y} \right)} = {{\delta {\ln \left( \frac{F_{0}}{F_{th}(N)} \right)}} - {2\frac{k^{2}\Delta x^{2}}{\omega_{0}^{2}}} - {2\frac{y^{2}}{\omega_{0}^{2}}\mspace{14mu} {or}}}},}} & (27) \\ {{{Z(y)} = {\frac{\Delta x\delta \sqrt{2}}{3\omega_{0}}\sqrt{{\ln \left( \frac{F_{0}}{F_{th}} \right)} - {2\frac{y^{2}}{w_{0}^{2}}}}\left\{ {{2{\frac{\omega_{0}^{2}}{\Delta \; x^{2}}\left\lbrack {{\ln \left( \frac{F_{0}}{F_{th}} \right)} - {2\frac{y^{2}}{w_{0}^{2}}}} \right\rbrack}} - 1} \right\}}},} & \left( 27^{\text{'}} \right) \end{matrix}$

where ΔX is the distance between two consecutive pulses, F0 is the peak fluence, F_(th) is the fluence of the ablation threshold and w0 is the radius of the point. X is given by v/PRR, where v is the sweep rate and PRR is the pulse repetition rate.

In cases where the focusing position of the beam is not corrected with increasing crater or line depth, which is usually the case, in the previous equation, ω₀ should be replaced by ω(z). The crater profile becomes a Z(y) “line” profile defined by:

Z(y)=Σ_(k=−N) ^(N) z _(k)(y).  (28)

To perform the analytical calculation of the previous equation, note that the number of pulses depends on y, i.e. it is not the same along the entire line profile. A number of pulses can be determined in the direction perpendicular to said line by the following equation:

$\begin{matrix} {{{N(y)} = {\frac{\omega_{0}}{\sqrt{2}\Delta x}\sqrt{{\ln \left( \frac{F_{0}}{F_{{th}{(N)}}} \right)} - {2\frac{y^{2}}{\omega_{0}^{2}}}}}}.} & (29) \end{matrix}$

It is important to note that the total number of pulses seen by a sample point for one pass is 2N+1. This expression leads to several comments. First of all, the line profile is no longer parabolic as it was the case for a plurality of static pulses. Second, depending on the value of Δ, the maximum number of pulses may be small, for example less than 10. In most materials, incubation effects must also be taken into account to obtain realistic results with the three previous equations.

When a material is ablated along a line, the fluence (F) expressed in J/cm2 can be expressed in the direction perpendicular to said line by the following equation:

$\begin{matrix} {{{F(y)} = {F_{0}{\exp \left( {- \frac{2k^{2}\Delta x^{2}}{\omega_{0}^{2}}} \right)}{\exp \left( {- \frac{2y^{2}}{\omega_{0}^{2}}} \right)}}},} & (30) \end{matrix}$

Δx describing the distance between the centres of two successive pulses, where F₀ is the peak fluence of the laser beam, k . . . , ω₀ is a minimum beam width,

Since the ablation threshold F_(th) is a parameter related to a fluence in J·cm⁻², the following assumptions about an ablation threshold value are preferred:

the redeposition of debris is negligible;

the ablation speed is not affected by the interaction with previous pulses;

there is no heat build-up or thermal effect on the surface of the material or plasma absorption;

there is no internal reflection.

With respect to the performance of a laser processing machine, a laser beam ablation rate can be defined by the following equation:

$\begin{matrix} {{\frac{dV}{df} = {\frac{f\delta \pi}{6}{{\ln \left( \frac{F_{0}}{F_{th}} \right)}\left\lbrack {{\ln \left( \frac{F_{0}}{F_{th}} \right)} - {\Delta x^{2}}} \right\rbrack}}}.} & (31) \end{matrix}$

A fluence ablation rate can also be defined by the following equation:

$\begin{matrix} {{\frac{dV}{df}/F_{0}} = {\frac{{fS}\; \Delta \; x}{F_{0}}.}} & (32) \end{matrix}$

A maximum ablation rate can then be defined by the following equation:

$\begin{matrix} {\left( {\frac{dV}{dt}/F_{0}} \right)_{m\alpha x} = {\frac{2}{e^{2}}{\frac{\delta}{F_{th}}.}}} & (33) \end{matrix}$

The equations defining the ablation rate have limitations because they do not allow the following parameters to be taken into account:

saturation due to the conicity of the machined part;

reflectivity (S,P);

incubation;

thermal effects.

The above parameters can be taken into account in a model for predicting the machining parameters for a target machining. By simulating the ablation profiles for all fluence ranges of the laser beam with the equations defining these ablation profiles and varying the parameters not taken into account in these equations. Such a method would be way too slow and would require too many experiments.

Perform an experimental design with a limited number of parameters and identify the parameters that have the greatest influence on the final results.

Determination of the threshold fluence using a diagonal scan model (Diagonal-Scan).

$\begin{matrix} {{\varphi_{th} = \frac{E_{p}}{e\; \pi \; \rho_{\max}^{2}}}.} & (34) \end{matrix}$

Consideration of the Local Angle of Incidence of the Beam

One way to introduce this local surface angle correction is to consider that the fluence varies with the distance to the crater centre r, since the angle between the beam direction and the surface varies from the crater rim to the crater centre. For a deep crater this angle is small on the walls but should reach 90° at the center. The conicity γ(r) of the walls is directly related to the derivative of the crater profile at the considered r:

$\begin{matrix} {{{\tan \left( {\gamma (r)} \right)} = {\frac{dz}{dr}{z_{k}(r)}}},} & (35) \end{matrix}$

which is equal to the slope of the tangent at each point of the crater profile (indicated in the figure below by the red line L). The local fluence after N pulses, F_N(r), can therefore be simply expressed as follows:

F _(N)(r)=F(r)cos(γ(r)).  (36)

It can be noted that the above equation can also be used even if the fluence remains constant from one pulse to the next, for example if the sample is already pierced or for a large crater shape.

Consideration of Non-Normal Beam Incidence

There are several applications that require cavity walls perpendicular to the surface, which is difficult to achieve if the laser beam is scanned normally towards the surface. However, negative or 0° wall conicities can be obtained if the laser beam is tilted in relation to the surface. To model this, the surface angle to the beam must also be considered. The previous equation thus becomes:

F _(abs,N)(r)=F _(N)(r)(1−R)cos(γ(r)+/β)  (37)

where β is the angle between the normal to the surface and the laser beam (FIG. 1). It should also be considered that the direction of ablation should be parallel to the direction of the beam rather than perpendicular to the original surface. The referential coordinates must therefore be adapted accordingly, as shown in FIG. 18 b.

As previously stated, function (27) or (27′) may be valid for groove ablation with the beam perpendicular to the surface, but will no longer be valid in cases where the beam is tilted with respect to the surface normal using processing parameters leading to negative conicities. In such cases, there will be a position y corresponding to more than one depth z. The minimum conicity that can be obtained with this approach therefore tends towards 0°.

To overcome this limitation, the model numerically performs this type of calculation by tilting the surface before calculating the depth (FIG. 17b ) rather than tilting the beam. The advantage of such a method is that it allows the simulation of machinings with negative conicity.

Consideration of Reflectivity

Polarization plays a major role in laser ablation processes of metals because of the role of the reflected beams. The evolution of reflectivity with the angle of incidence α is given by the Fresnel formulae for s, p circular polarizations:

$\begin{matrix} {R_{s} = {\frac{\sqrt{n^{2} - {\sin^{2}\alpha}} - {\cos \; \alpha}}{1 - n^{2}}}^{2}} & (38) \\ {R_{p} = \left| \frac{{n^{2}\cos \; \alpha} - \sqrt{n^{2} - {\sin^{2}\alpha}}}{{n^{2}\cos \; \alpha} + \sqrt{n^{2} - {\sin^{2}\alpha}}} \right|^{2}} & (39) \\ {R_{circul\alpha r} = \frac{R_{s} + R_{p}}{2}} & (40) \end{matrix}$

The fluence absorbed by the material is then given by the following expression:

F _(abs)(r)=F(r)(1−R)  (41)

where R is R_(S), R_(P) or R_(circular) depending on the direction of polarization radiation. If we introduce into the correction described by the equation allowing the contribution of the surface reflectivity to be taken into account, the local fluence absorbed after N pulses then becomes:

F _(abs,N)(r)=F _(N)(r)(1−R)cos(γ(r))  (42)

In the case of the referential, we use (FIG. 18b ), α=90°−γ. In general, for metals, for low conicity, i.e. for a high value of a, the p-polarized reflectivity becomes low and the transmission is improved, so the ablation efficiency is also improved for a low conicity. For the oblique beam, the same calculation can be done with α-β instead of β. With a medium polarization, the effect is taken into account with a multiplicative coefficient at F₀(0) in the equation defining the local fluence after N pulses using the reflectivity equations of R_(s) and R_(p). This coefficient is 1 for conicity greater than 20°, and can reach 1 for conicity less than 20° also depending on the value of β.

Compared to what is proposed in the prior art, the simulation method of the invention comprises that each material whose material laser interaction parameters have been modeled and added to the database is analyzed by ellipsometry to determine its complex refractive index (n+ik) as a function of the wavelength of the laser radiation.

Thus, the reflectivity is taken into account as follows:

the reflectivity is calculated as Rs, Rp or (Rs+Rp)/2 according to the Fresnel formulas of equations (38), (39), and (40), for S, P or circular polarization, respectively;

the resulting R, which will be between 0 and 1, is subtracted from the incident fluence, giving an absorbed fluence given by:

F _(absorbed)(γ)=F _(c)(0)×sin(γ_(n)(γ)+β)×(1−R).

Gaussian Beam Propagation

Since laser treatments are usually performed with the (Gaussian) beam focused using a lens with a certain focal length, the beam propagation must be taken into account as most experiments are performed without adjusting the surface distance to the focal plane. The depth of ablation increases. The variation of the beam radius as a function of the focal plane distance is given by:

$\begin{matrix} {{\omega (z)} = {\omega_{0}\sqrt{1 + \left( \frac{z}{z_{R}} \right)^{2}}}} & (43) \end{matrix}$

where z_(R), is the Rayleigh distance, given by the relation:

$\begin{matrix} {z_{R} = \frac{\pi \omega_{0}^{2}}{\lambda}} & (44) \end{matrix}$

In the case of a non-normal incidence, the fluence distribution as a function of z and the radial coordinate r can be determined numerically by applying a geometric correction as shown in FIG. 19. The centre of the beam with an angle of inclination α=90°−β can be described as follows:

z _(R) =r tan(β)+f,  (45)

where f is the distance from the surface to the focal plane. On the other hand, the line perpendicular to this central axis of the beam and passing through the generic point P (r, z) is given by:

z _(P) =r tan(β+90°)+O _(p)  (46)

To determine the fluence at the generic point P(r,z), it is necessary to know the Gaussian distribution of the intensity along the plane z_(P) and, in turn, the distance between the point BP(r,z) and the focus f. The point BP can be calculated from the intersection of the two lines z_(B) and z_(P):

z _(B) =z _(P)  (47)

Z _(BP) =r _(BP) tan(β+90°)+O _(P)  (48)

O _(P) =z _(BP) −r _(BP) tan(β+90°)  (49)

The distances z′ and r′ can then be easily obtained from the Cartesian coordinates of the points BP and P:

z′=√{square root over ((r _(BP))²+(z _(BP) −f)²)}  (50)

r′=√{square root over ((r _(BP) −r _(P))²+(z _(BP) −z _(P))²)}  (51)

From the two distances z′ and r′, it is then possible to calculate the local fluence, using the width of the beam on this plane described by the line z_(P):

$\begin{matrix} {{{F\left( z_{P} \right)} = {F_{0}{\exp \left( {2\frac{{r^{\prime}}^{2}}{\omega \left( z^{\prime} \right)}} \right)}}},{{o\; ú\; {\omega \left( z^{\prime} \right)}} = {\omega_{0}{\sqrt{1 + \left( \frac{z^{\prime}}{z_{R}} \right)^{2}}.}}}} & (52) \end{matrix}$

Since laser treatments are usually performed with the (Gaussian) beam focused with a lens of a certain focal length, the beam propagation must be taken into account, as most experiments are performed without adjusting the surface distance to the focal plane as the depth of ablation increases.

Contrary to the prior art, the model included in the algorithm implemented by the central unit takes into account the beam propagation. The variation of the beam radius as a function of the focal plane distance is given by:

$\begin{matrix} {{\omega (z)} = {\omega_{0}\sqrt{1 + \left( \frac{Z}{z_{R}} \right)^{2}}}} & (53) \end{matrix}$

where z_R, Rayleigh's distance, is:

$\begin{matrix} {z_{R} = \frac{\pi \omega_{0}^{2}}{\lambda}} & (54) \end{matrix}$

For example, this results in the following beam distribution profile in the focal plane area for a beam with a focal diameter of 30 μm and a wavelength of 1030 nm. In this case, the Rayleigh distance is 686.27 μm.

Simulation of Machining Using a Precession Laser Beam

Precession laser ablation is a technique that allows to perform straight-sided cutting and drilling at high-speed that is more advantageous and efficiently than, for example, using a nozzle and a gas or scanning the beam at high speed with a scanner head. The precession movement can be achieved by using a rotating mirror or other rotating element in the optical path of the laser. With a rotating mirror, it is easy to control the rotational speed of the beam co (in rpm), the positions of the center of rotation of the beam BFG (best focus-global beam) and the focal plane of the beam BFI (best focus-individual beams) relative to the surface and the angle of attack B.

In this ablation configuration, the pulse superposition will be different from the ablation produced by a linearly scanned beam. To determine the ablation profile of a drilled hole, the distance between spots dx, which depends on the precession radius r_(P), and the angle θ, must be known.

The precession radius r_(P) is given by:

$\begin{matrix} {r_{P} = {\left( {\frac{BFG}{\cos \beta} - z} \right)\tan \beta}} & (55) \end{matrix}$

where BFG is the distance to the centre of rotation of the beam, β is the angle of attack and z is the distance to the material surface.

The distance between two consecutive pulses dx can be calculated as the distance between two vertices of a polygon with NR sides (eq. 57), where NR is the total number of pulses during one revolution. NR depends on the pulse repetition rate PRR (Hz) and the rotational speed co (rpm) (eq. 56).

$\begin{matrix} {{NR} = \frac{PRR}{{\omega/6}0}} & (56) \\ {{dx} = {2r_{P}{\sin \left( \frac{n\theta}{2} \right)}}} & (57) \end{matrix}$

Thus, there is a superposition between pulses if the distance dx between the first pulse no and the pulse n_(N) is less than the radius of an ablated crater r_(c) (FIG. 20a ), which is given by:

$\begin{matrix} {r_{c} = {w_{0}\sqrt{\ln \left( \frac{F_{0}}{F_{th}} \right)}}} & (58) \end{matrix}$

Finally, the depth of ablation produced by each pulse n during one revolution as a function of the distance to the center is given by:

$\begin{matrix} {z_{n} = {- {\delta \left\lbrack {{\ln \left( \frac{F_{0}}{F_{th}\left( {n + {N{R/2}}} \right)} \right)} - {2\frac{n^{2}dx^{2}}{w^{2}}} - {2\left( \frac{\left( {y + y_{C}} \right)^{2}}{w^{2}} \right)^{p}}} \right\rbrack}}} & (60) \end{matrix}$

where δ is the depth of penetration of the radiation into the material, F₀ is the maximum fluence, F_(th) is the ablation threshold, w is the radius of the Gaussian beam and p is the Gaussian order. This depth must be iterated R times, where R is the total number of revolutions, given by:

R=proc.time*ω/60.  (61)

where proc.time is the total machining time. Note also that the radius of the Gaussian beam depends on the distance to the focal point BFI. In the equations described above, it is important to note the following: the w will change as the machining depth increases. After each ipulse the r_(P) will change if part of the machining has already been ablated.

Preferably, the formulas described above can be used individually or in combination to implement an algorithm for simulating a machining profile in two and or three dimensions.

BRIEF DESCRIPTION OF THE FIGURES

These and other aspects of the invention shall be clarified in the detailed description of particular embodiments of the invention, reference being made to the drawings of the figures, in which:

FIG. 1 represents the main characteristics of a beam having a Gaussian intensity profile;

FIG. 2 represents a matrix of craters with increasing numbers of pulses and pulse energies;

the FIG. 3a -FIG. 4b represent: FIG. 3a and FIG. 4a , a plurality of pulses along a line;

FIG. 3b , a fluence profile obtained transversely to the line described in FIG. 3a ; FIG. 3c represents a distribution of pulse number N transversely to the line described in FIG. 3a ; FIG. 4b , the groove depth obtained transversely to the line described in FIG. 4 a;

FIG. 5 shows simulation results obtained fora cross section of a line representing an ablation rate (left) and a maximum ablation rate as a function of the fluence of the laser;

FIGS. 6a and 6b show in FIG. 6a the displacement of the sample in the direction of the axes Y and Z for the implementation of the D-Scan method and in FIG. 6b the ablation lobes resulting from the implementation of the experimental protocol of FIG. 6 a;

FIGS. 7a and 7b show a schematic representation of the experiments to be performed to estimate the value of δ using the skin depth method;

FIGS. 7c and 7d show the experimental results compared to the models for the experiments shown in FIGS. 7a and 7b respectively;

FIGS. 8a and 8b show the results for the dimensions ρ_(min) and ρ_(max) from the experiments shown in FIGS. 7a and 7b for laser beam pulse frequencies of 5 KHz and 100 KHz respectively.

FIG. 9 shows the results of the 5 KHz and 100 KHz experiments of FIGS. 8a and 8b in a graph showing the fluence F_(th) as a function of pulse number N;

FIG. 10 shows a simulation of the measured ablation as a function of the angle of incidence for a polarization Son the left and P on the right;

FIG. 11 shows the propagation of a Gaussian beam and represents the beam diameter w as a function of a propagation distance z;

FIGS. 12 and 13 represent experimental and simulation results for two values of δ. FIG. 12 represents the maximum depth as a function of the beam fluence. FIG. 13 represents the maximum depth as a function of scan speed of the beam;

FIGS. 14a, 14b, 15a, 15b, and 15c represent experimental and simulated profiles obtained on stainless steel 316L;

FIGS. 16a and 16b show in FIG. 16a a cross-section of a scan and the corresponding simulated cross section (FIG. 16b ) obtained on stainless steel 316L;

FIGS. 17a and 17b show a schematic representation of a machining with a normal angle of incidence in FIG. 17a and a machining with a non-normal angle of incidence in FIG. 17 b;

FIG. 18 shows a schematic representation of the angles of a light beam in experiments for the calculation of the local fluence at a non-normal incidence;

FIGS. 19a and 19b show a comparative example between a simulation of a machining profile according to the prior art (FIG. 19a ) and a simulation according to a particular embodiment of the method of the invention for a machining carried out with a machining laser beam having a non-zero angle of attack with the normal to the surface of the material to be machined.

FIG. 20a shows a schematic representation of machining with a precessing laser machining beam;

FIGS. 20b, 20c, 21b and 21d show simulations of machining according to a particular embodiment of the method of the invention;

FIGS. 21a and 21c show machining results in the form of cross sectional images of a machining, the machining simulations illustrated in FIG. 21b and FIG. 21c are superimposed to the machining results respectively;

FIG. 22a shows a schematic representation of the points of impact of a precessing or rotating machining laser beam, FIG. 22b shows a detail of FIG. 22 a;

FIGS. 23 to 31 show flow diagrams of preferred embodiments of the invention.

The drawings of the figures are not to scale. Generally, similar elements are denoted by similar references in the figures. The presence of reference numbers in the drawings may not be considered as limiting, even when such numbers are indicated in the claims.

EXAMPLE OF AN EXPERIMENTAL DEVICE FOR IMPLEMENTING THE METHOD OF THE INVENTION

The experiments carried out to compare the models described above were performed on samples of polished stainless steel 316L and 316 and TiCr6Sn4. The tests were performed in air using a Satsuma HP2 (Amplitude Systems) femtosecond laser with a pulse duration of about 330 fs, a radiation wavelength of 1030 and a maximum power of 20 W at 500 kHz. The beam was focused on the surface of the samples using a telecentric lens with a focal length of 100 mm, producing a spot radius of about 10 μm determined using the D² method.

Morphological and topographical analysis of the processed samples was performed using a confocal optical microscope (Olympus LEXT OLS4100).

The D² method was used to calculate the threshold fluence (F_(th)) values. To determine the value of δ, several line scans with increasing pulse energy were produced and their depths measured. The δ value applied in the model was varied until the best match with the experimental results was obtained.

Preferably, the greater the amount of results provided to the database, the more accurately the laser machining parameters are determined.

DETAILED DESCRIPTION OF CERTAIN EMBODIMENTS OF THE INVENTION

FIG. 1 shows an intensity profile of a Gaussian laser beam. The center of the beam is defined by the value 0 along the abscissa axis defining a beam radius r and a maximum intensity I₀ is observed at the center of the beam. A beam radius coo relative to the beam center is defined for a beam intensity value equal to I₀/2.

FIG. 2 shows a crater matrix to be machined in order to implement the ablation threshold determination by the D² method. FIG. 2 describes a preferred embodiment for producing a crater matrix with the laser machining machine of the invention. The matrix shown in FIG. 2 shows a number of pulses for each of the columns between 10 and 110 and beam energies between 10% and 100% of a maximum beam value that can be provided by the laser machining machine of the invention. The realization of such a matrix of craters while respecting the beam energies and the numbers of pulses indicated, followed by the topological and/or dimensional characterization of the craters of this matrix thus enables to define an ablation threshold, an incubation coefficient or even a beam radius. Preferably the production of the crater matrix should be performed with increasing pulse numbers and pulse energies, as illustrated by the arrows in FIG. 2. In order to have sufficient experimental points, at least 20 energy values and several pulse number values less than 100 pulses are required. The diameter of each crater can be measured from images obtained with an optical or electron microscope or a profilometer.

FIG. 3a shows a plurality of pulses along a line. A pulse is represented by a circle with its center on a line. The centres of each circle, i.e. the points where the maximum intensity of each pulse is maximum I₀, are spaced along a line of: Δx=v/PRR. v being the speed of movement of the beam and PRR being a number of pulses per unit time. FIG. 3a shows that for pulses along a line, with a distance between each circle less than the radius of each circle, then statistically more pulses reach the centre of the line than the edges of the line.

FIG. 3b shows a fluence profile obtained transverse to the line described in FIG. 3a . This fluence profile shows that the fluence at the center of the line is much higher than at the edges of the line.

FIG. 3c shows a distribution of the pulse number N transversaly to the line described in FIG. 3a . As expected from the description in FIG. 3a , statistically, a greater number of pulses reach the center of the line than the edges of the line.

FIG. 4a shows as in FIG. 3a a plurality of pulses along a line. FIG. 4b , shows the groove depth obtained transversaly to the line described in FIG. 4a with the following experimental parameters:

material: steel;

ω₀=12.5 μm;

F=3 J/cm²;

V=50 mm/s;

PRR=200 kHz.

FIGS. 7c and 7d represent the experimental results compared to the corresponding models obtained by the experiments described in FIGS. 7a and 7b . FIG. 7c shows the maximum scan depths obtained on a steel substrate as a function of the maximum fluence. FIG. 7d shows the maximum depths of scans obtained on a steel substrate as a function of the offset speed or the scan speed.

FIGS. 8a and 8b show point clouds obtained for experiments described in FIGS. 6a and 6b for which a regression according to an algorithm of the invention allows to determine the parameter ω₀.

FIG. 9 shows the experimental results of the D-scan method for two pulse frequencies of 5 KHz and 100 KHz. The regressions for these two sets of experimental results are also illustrated by the lines. The parameters of these two regressions are indicated by Gama, Beta and Alfa.

FIG. 10 shows a parameter for the interaction of a polarized light beam with a steel substrate having the following properties: n=0.9+2.25i. FIG. 10 shows the interraction parameter for incidence angles between 0° and 90°, for a polarization s and for a polarization p.

The curve representing the polarization p is the one with a maximum at 100% for an incidence angle of 90°.

FIG. 11 shows a simulation of the beam width w_(z) at its minimum width. The minimum beam width is shown at a distance z in a propagation direction z having a value of 0. The simulated beam profile shows the divergence of the beam in the direction z.

FIG. 12 shows experimental results and simulations for two values from δ concerning maximum depth as a function of beam fluence. For peak fluence values less than 2 J/cm² and greater than 4.5 J/cm², the experimental results vary from the two simulations having a value δ of 32 nm and 34 nm. This figure illustrates very well the need to have a simulation algorithm that is sufficiently complete to take into account as many machining parameters as possible. Thus, a good taking into account of the material laser interaction parameters allows to obtain a good correlation between a machining result and the predictions of a simulation for a constant delta.

FIG. 13, showing the maximum depth as a function of beam scanning speed, shows a good correlation between the experimental ablation results and the beam scanning speed. The experimental results also show a relatively good robustness when varying the parameter of δ from 32 nm to 34 nm.

FIGS. 14a and 14b show experimental and simulated profiles obtained on stainless steel 316L with a maximum fluence of 2.14 J/cm², sweep rate 230 mm/s (FIG. 14a ) and 1.25 J/cm² and sweep rate 50 mm/s (FIG. 14b , p-polarization and 250 kHz frequency). A good correspondence between experimental and simulated results is observed with respect to the maximum depth of ablation as well as the ablation profile. On the other hand, the redeposition of material at the edges of the ablation is not taken into account by the simulation. The simulation method of the invention could allow taking into account the redeposition of material on the edges of the ablation experimentally performed by developing an adapted algorithm.

FIGS. 15a, 15b and 15c represent experimental and simulated profiles obtained on stainless steel 316L with a maximum fluence of 1.92 J/cm², scanning speed 100 mm/s, frequency 100 kHz, the FIG. 15a is obtained after 17 scans, the FIG. 15b after 33 scans and the FIG. 15c after 65 scans. The same observations as for FIGS. 14a and 14b can be observed.

FIGS. 16a and 16b show in FIG. 16a a cross section of a scan and the corresponding simulated scan (FIG. 16b ) obtained on stainless steel 316L with a maximum fluence of 2.9 J/cm², scanning speed of 500 mm/s, frequency of 100 kHz and 500 passes with the beam tilted of 15° in relation to the normal to the surface. The model of the invention allows a good taking into account of the angle of incidence of the beam incident on the material to be machined. As shown in FIG. 16a , there is a good correspondence between the micrography of the machining (ablation) and the simulation which is superimposed.

FIGS. 17a and 17b show machining profiles obtained with an angle of incidence normal to the surface of the workpiece and with an angle of incidence describing an angle 90°-β in relation to the surface of the workpiece. For example, these schematic representations with different angles of incidence in relation to the surface of the workpiece enable to calculate the local fluence.

FIG. 18 shows a schematic representation of the angles of a light beam in experiments for the calculation of the local fluence under a non-normal incidence. In particular, this schematic representation illustrates the divergence of the beam having a non-normal incidence with a substrate to be machined. The peculiarities caused by this divergence of the focused beam are thus illustrated.

FIGS. 19a and 19b show a comparative example between a simulation of a machining profile according to the prior art (FIG. 19a ) and a simulation according to a particular embodiment of the invention (FIG. 19b ) for a machining performed with a machining laser beam having a non-zero angle of attack with the normal to the surface of the material to be machined. Both FIGS. 19a and 19b are multipass profiles obtained on stainless steel (in both cases δ=32 nm, F_(th)(1)=0.1 J/cm², S=0.8, n=0.9-2.25i) for w₀=14 μm, PRR=200 kHz, Δx=0.75 μm, F₀=8 J/cm² and an angle of incidence of 8°. The numbers in the figures show the number of passes of the machining laser beam for each of the simulated machining profiles in two-dimensions. The simulations were therefore carried out for 50, 100, 200, 300, 402 passes. The particular embodiment of the invention which allows the processing of machining laser beams having a non-zero angle in relation to the surface of the material allows a very representative simulation of how the machining would be using the same laser parameters. In the case of the simulation using prior art methods, it is observed that a high number allows a machining flank to be approached perpendicular to the surface of the material but without ever reaching it, which is problematic for the simulation of machining with a conicity of the flanks in relation to the surface of the material which is well controlled. Indeed, the precession or other machining laser beam modules that allow machining with a well-controlled beam angle require suitable simulation tools. The particular embodiment of the invention seems to be one of them. The ordinate axis in FIGS. 19a and 19b represents a simulated ablation depth expressed in μm. The abscissa axis of FIGS. 19a and 19b represents a dimension y corresponding to a dimension parallel to the surface of the material to be machined.

FIG. 20a shows a schematic representation of a machining with a precessional laser machining beam. It shows a material to be machined whose surface corresponds to the coordinate z=0 in relation to the axis z essentially perpendicular to the surface of the material. The laser beam is rotating, especially during precession as indicated by the circular arrow indicating a direction of rotation. The beam comes from the side of this circular arrow and is directed towards the material to be machined. The machining laser beam is shown as passing through the material so as to illustrate the trajectory of the precession laser beam. Thus the laser beam describes a point BFG located on a plane BFG essentially parallel to the surface of the material. As can be seen, the precession beam always passes through the same point at the point BFG. Furthermore, the plane BFI is the plane in which the precessional laser beam is focused. The precession laser beam rotates around an axis represented as being normal to the surface of the material. The precession laser beam therefore describes an angle β for all precession positions of the beam with respect to the axis of rotation of the beam.

FIGS. 20b, 20c, 21b and 21d show machining simulations according to a particular embodiment of the invention for a precession laser beam as shown schematically in FIG. 20a . FIGS. 20b, 20c, 21b and 21d show machining profiles in two-dimensions along the axes y, z, with the ordinate axis representing the dimension z (depth of ablation) and the abscissa axis representing the dimension y. FIG. 20b corresponds to a simulation according to the particular embodiment of the invention allowing the simulation of a precession laser machining beam with BFG=150 μm and BFI=−50 μm. Similarly, FIG. 20c was obtained with BFG=−100 μm and BFI=−500 μm. The simulations in FIGS. 20b and 20c were obtained with the laser parameters: w₀=10 μm, PRR=500 kHz, E_(p)=20 μJ, ω=30 000 rpm, 13=4.5°, machining time=50 ms. FIG. 21b corresponds to a simulation according to the particular embodiment of the invention allowing the simulation of a precession laser machining beam with BFG=0 μm and BFI=−400 μm. Similarly, FIG. 21d was obtained with BFG=−100 μm and BFI=−500 μm. The simulations in FIGS. 21b and 21d were obtained with the laser parameters: w₀=12 μm, PRR=100 kHz, E_(p)=27 μJ, ω=30 000 rpm, β=4.5°, machining time=50 ms.

FIGS. 21a and 21c show machining results in the form of cross-sectional cutting images of a machining, the machining simulations shown in FIG. 21b and FIG. 21c are superimposed on the machining results respectively. It is thus observed that the machining simulations allow a good anticipation of the machining obtained on the basis of the same parameters for a precession laser machining beam.

FIG. 22a shows a schematic representation of the points of impact of a precessing or rotating machining laser beam, FIG. 22b shows a detail of FIG. 22a . FIG. 22a , shows the pulse superposition of a precession laser machining beam. Around an axis represented here by a point because it is a projection of the laser pulses in a plane, preferably on the surface of the material. The axis of rotation of the precession laser beam is therefore perpendicular to the surface of the material and therefore represented by a point. Each circle represents one pulse. The center of each pulse is located on a precession circle centered on the axis of rotation of the machining laser beam. The precession circle has a radius r_(p). The distance between the center of each adjacent circle on the precession circle is dx. Preferably dx is constant during machining. The center of each adjacent circle describes an angle θ in relation to the precession axis. FIG. 22b represents a detailed view of FIG. 22a in which, the center of the circles corresponding to a pulse is named n₀, n₁, n₂, n₃, . . . . The distance between n₀ and n₁ is d₁, the distance between n₀ and n₂ is d₂, the distance between n₀ and n₃ is d₃.

FIGS. 23 to 31 represent flow charts describing the various elements of a laser machining system and the possible steps or decisions that can be taken. In particular, these flowcharts show the means by which information can be channeled and lead to the machining of a material with optimum parameters or to the enrichment of a material database and/or functions learned by machine learning.

FIG. 23 represents a Functional Graph for Step and Transition Control (Grafcet) or functional flowchart of the different actions available according to the invention. In particular, this flowchart shows that the method of the invention may decide, in the event that the parameters of a material are not available in a database, to resort to the launching of a predefined test machining in order to allow the acquisition of data in relation to said material to be machined.

FIG. 23 shows the following procedure: A material to be machined is communicated to the laser machining system, which queries whether information about this material is present in a material database or in a database of learned functions. In the case of a positive answer, it is then possible to choose between manual or automated determination of the machining parameters.

If the machining parameters are determined in an automated manner, the learning system sends a request for knowing the machining result sought. This result to be achieved can also be communicated when the material to be machined is communicated. On the basis of the result to be achieved (target result) defined by the user, the laser machining system and in particular the central unit collects information about the characteristics of the machining device as well as the information available in one of the databases about the information of the material to be machined. The learning central unit then enables the exploitation of the means for modelling and/or machine learning by taking into account the information of the material to be machined, the characteristics of the machining device and the machining result sought. The learning central unit can then generate optimal machining parameters. These optimal machining parameters are then transmitted to the laser machining device in order to start the laser machining of the material to be machined according to the result to be achieved. It is possible that the learning central unit is connected to the laser machining device via a network connection so that the learning central unit can be relocated in relation to the laser machining device.

If the machining parameters are determined manually, the central unit sends a request to obtain parameters from the user to obtain a result to be achieved. Preferably the parameters transmitted by the user are those parameters which the user believes to be the optimum results. However, several iterations are often necessary for the operator of the laser machining system to specify parameters allowing to obtain the result to be achieved. When parameters are communicated by the user, then modeling means are implemented. The modelling means preferably comprises a model including an algorithm allowing an estimation of the machining on the basis of the parameters provided. Thus simulation means coupled to the modeling means enable to carry out a simulation of the result expected to be obtained on the basis of the parameters communicated. The operator can then compare the simulation of the result on the basis of the parameter with the result to be achieved. If the operator considers that the simulation does not conform (or sufficiently conforms) to the result to ba achieved and communicates this to the central unit, then the central unit offers a choice between manual or automated determination of the machining parameters. The operator can then decide to test several different machining parameters manually until a simulation of the result expected on the basis of the communicated parameters conforms to the result to be achieved. When the operator is satisfied with the simulation of the result (machining) expected on the basis of the communicated parameters, then he can decide to start machining on the basis of the parameters used in the last modeling and simulation. The operator can choose at any time to use the automated parameter search so that the laser machining system uses the automated mode as described above. The central unit then sends a request to know the machining result sought and determines the optimal machining parameters as described above and as represented by the flowchart. During the modeling step in manual mode, the modeling means have access to the material database as well as to the characteristics of the machining device.

If the query as to whether information on the material to be machined is present in a material database or in a database of learned functions is answered in the negative, the laser machining system then defines predefined machining tests. These predefined machining tests enable the laser machining system to produce a predefined machining on a material, preferably the material in the material database being identical to the material to be machined, the machining being analysed by the analysis unit and the analysis results being transmitted to the central unit.

The results of the analysis are either communicated to modeling means for extracting physical or material light interaction parameters specific to the material to be machined. Modeling means for receiving information about the characteristics of the machining device. The parameters determined by the modelling means are then machining parameters which are communicated to the material database in order to enrich it. When the material to be machined is known in the database, it is then possible to continue the machining process and in particular towards the stage proposing a manual (this invention) or automated parameter search.

When the analysis results are reported to the central unit, the central unit of the machining system having access to the features of the machining system allows information about the material being machined to be generated and stored or updated in the material database.

The preferred embodiment described in the previous paragraph is shown in FIG. 24. Thus, the same database can be used to group together, for the same material, the physical characteristics of the materials as well as the information relating to the laser machining system so as to form a learning database for the implementation of an machine learning algorithm. The flow chart described for FIG. 23 applies to FIG. 24 mutatis mutandis.

FIG. 25 shows an embodiment of the invention for determining the optimum machining parameters in a deterministic regime for a laser machining system. In particular, this embodiment is based primarily on the use of known algorithms as described by the invention. This embodiment also allows the simulation of ablations which would be obtained for laser parameters obtained by the algorithms of the invention or defined by a user. This figure describes the embodiments detailed in FIGS. 26 to 28. These embodiments are derived from the flowcharts shown in FIGS. 23 and 24.

FIG. 26 shows an embodiment of the invention allowing the feeding of physical material characteristics to the material database, for example in the case where a predefined test has been machined and for which the results are communicated to the modeling means as described for FIGS. 23 and 24. The arrow connecting the material database to the modelling means can be interpreted in either direction. The purpose of this embodiment is the determination and the storage of experimental parameters in a database. From a machining performed by a laser machining device, the results of this machining are communicated to a central unit which allows the experimental data to be processed in order to extract physical parameters of materials.

FIG. 27 shows an embodiment of the invention allowing the simulation without necessarily requiring the use of a laser machining device. Physical characteristics of materials previously stored in the material database allow a simulation of the results that could be obtained for machining according to user parameters and characteristics of the machining device.

FIG. 28 shows an embodiment of the invention for determining optimal machining parameters from the results to be achieved, the characteristics of the laser machining device and from a database comprising physical characteristics of materials, for example observed and stored in the material database during previous machinings.

FIG. 29 shows an embodiment of the invention and in particular an embodiment of the invention based on machine learning means for taking into account physical characteristics of materials or functions learned by the machine learning. A learned function is the result of the learning of a learning algorithm itself included in a learning function. This information is preferably stored in a database or in several databases. This embodiment is preferably based on a non-deterministic optimal parameter determination regime, i.e. not necessarily using modelling means comprising an algorithm defining deterministic functions.

Preferably, the embodiment detailed in FIG. 29 includes a learning function whose learned functions are stored in a database. This embodiment allows the method of the invention to learn and refine the optimal parameters based on past experiences and observations. For example, experiments carried out on other laser machining devices are taken into account in learning the algorithm for determining optimal machining parameters. This embodiment also allows for a specific experiment to be carried out to collect experimental data in order to obtain optimal machining parameters for a target machining. For example, predefined tests are available in the laser machining system of the invention and are applicable when it is deemed necessary to obtain machining with high fidelity with respect to a target machining.

FIG. 30 shows a special implementation of FIG. 29 allowing automated learning and recording of learning functions as learned fonctions. FIG. 30 shows the feeding of the database of learned functions by the machine learning in non-deterministic mode. Preferably, this embodiment implements a predefined machining test, which can be started on a selected material. For example, a predefined machining test is defined in FIGS. 2, 3 a, 6 a, 6 b, 7 a, 7 b. These tests allow, after characterization of the machinings obtained by the analysis unit, the machine learning and the determination of a learned function. The results of these tests and the learned functions generated by taking these tests into account are then stored in a database.

FIG. 31 shows an embodiment for machining a workpiece using the optimum parameters generated by the central unit or by the learning means on the basis of the learned functions corresponding to the material to be machined. This embodiment allows to obtain machining results very close to the desired target machining and very quickly. The FIG. 31 shows the prediction of the optimum machining parameters in non-deterministic regime.

The present invention has been described in relation to specific embodiments, which have a purely illustrative value and should not be considered as limiting. In general, the present invention is not limited to the examples illustrated and/or described above. The use of the verbs “comprise”, “include”, “consist”, or any other variant, as well as their conjugations, can in no way exclude the presence of elements other than those mentioned. The use of the indefinite article “a”, “an”, or the definite article “the” to introduce an element does not exclude the presence of a plurality of such elements. Reference numbers in claims shall not limit their scope.

In summary, the invention can also be described as follows.

A method for simulating laser machining of a material by a laser machining system comprising the following steps:

a) providing a central unit with:

-   -   information relating to the material to be machined: delta δ,         threshold fluence, incubation coefficient S, complex refractive         index n+ik,     -   information about the laser machining system:         -   information relating to polarization,         -   a pulse energy E_(P),         -   a diameter of said machining laser beam at a focal point w,         -   an order from a Gaussian p,         -   a pulse repetition rate PRR n,         -   a wavelength;

b) determining with said central unit on the basis of:

-   -   said information relating to said material to be machined, and,     -   to the laser processing system,

a machining profile in two-dimensions corresponding to the simulation of a machining of said material to be machined with said laser machining system. 

1. A method for simulating a laser machining of a material by a laser machining system capable of emitting a machining laser beam, said method comprising the following steps: a) providing a central unit with: information relating to the material to be machined comprising: a delta δ relative to a penetration depth of the machining laser beam in the material, a threshold fluence F_(th) relative to a minimum energy density of the laser beam allowing a machining, an incubation coefficient S, between 0 and 1, relating to a dependence of the threshold fluence F_(th) on a number of pulses of the laser beam at a single point, a complex refractive index n+ik, information relating to said laser machining system comprising: information relating to a polarization of said machining laser beam, a pulse energy of said machining laser beam E_(P), a diameter of said machining laser beam at a focal point w, an order of a Gaussian p of said machining laser beam between 1 and 20, a pulse repetition rate PRR of pulses n of the machining laser beam, a wavelength of the machining laser beam; b) determining with said central unit on the basis of: said information relating to said material to be machined, and, said information relating to said laser machining system, a machining profile in two and/or three dimensions corresponding to the simulation of a machining of said material to be machined with said laser machining system.
 2. The method according to claim 1 wherein: the laser machining system is a system capable of emitting a laser machining beam moving in a direction x substantially parallel to a machining surface of said material, and in that, the step (a) further comprises providing the following information relating to said laser machining system: a distance of the focus point of the laser machining beam from the surface of the material to be machined, a speed (v) of movement of the laser machining beam relative to said material to be machined, an angle of incidence of the machining laser beam with respect to said surface of the material to be machined, a number of lines to be machined, a distance between said lines previously defined, number of passes of the machining laser beam on each line to be machined.
 3. The method according to claim 1 wherein: the laser machining system is a system capable of emitting a laser machining beam in rotation about an axis of rotation substantially perpendicular to a machining surface of said material, and in that, the step (a) further comprises providing information on said laser machining system: a rotational speed ω of the laser machining beam, a distance BFG from a surface of said material to be machined for which the rotating laser beam describes a fixed spot for all rotating laser machining beam positions, a distance BFI from a surface of said material to be machined at which the rotating laser beam is focused, an angle of incidence β of the rotating machining laser beam for all positions of the rotating laser beam with respect to the normal to the surface of the material, an activation time of the rotating machining laser beam.
 4. The method according to claim 3 wherein the machining laser beam is a machining laser beam in rotation about an axis of rotation and a point of rotation located at said distance BFG from said surface of said material to be machined.
 5. The method according to claim 3 wherein step b) of determining said machining profile in two and/or three dimensions comprises a step of determining a precession radius r_(p) of a machining of said material by said rotating machining laser beam by the following formula: $r_{P} = {\left( {\frac{BFG}{\cos \beta} - z} \right)\tan {\beta.}}$
 6. The method according to claim 5 wherein two successive pulses n of the rotating machining beam are separated by a distance dx, and wherein step b) of determining said machining profile in two dimensions comprises a step of calculating an ablation depth a by the following formula.
 7. The method according to claim 1 wherein step b) of determining said machining profile in two dimensions comprises a step of determining an ablated crater radius r_(c) for a pulse of the machining laser beam by the following formula: $r_{c} = {w_{0}{\sqrt{\ln \left( \frac{F_{0}}{F_{th}} \right)}.}}$
 8. (canceled)
 9. The method according to claim 1 wherein said machining profile is defined in two dimensions according to: an axis y representing a direction substantially parallel to the surface of the material to be machined, and, an axis z representing a direction substantially perpendicular to the surface of the material to be machined, the axis z corresponding to an ablation depth z_(n) with respect to the surface of the material.
 10. The method according to claim 1 wherein delta δ is a constant parameter.
 11. The method according to claim 1 wherein said central unit is configured to execute a computer program for the determination of said machining profile in two and/or three dimensions.
 12. The method according to claim 11 wherein said central unit configured to execute said computer program is accessible on dematerialized storage and computing means, for example a cloud, in particular a platform as a service and more particularly a software as a service.
 13. The method according to claim 1 wherein said information relating to said material to be machined is derived from a material database, wherein each material comprises the following information: a delta δ, a threshold fluence Fth, an incubation coefficient S, a complex refractive index n+ik, and wherein said central unit comprises communication means for communicating with said material database.
 14. The method according to claim 13 wherein said central unit comprises simulation means so as to be able to determine said machining profile in two dimensions on the basis of information relating to said material to be machined from said database, and on the basis of information relating to said laser machining system.
 15. The method according to claim 1 wherein the method further comprises the following additional steps: c) providing a laser machining system comprising: said material to be machined; said laser machining system comprising: a laser source for emitting a laser machining beam onto said material to be machined; a control unit for controlling the emission of said machining laser beam from said laser source; said central unit for controlling said control unit; said communication means for communicating with said database; communication means for communication from said central unit to said database; d) machining said material with said laser source configured with said laser machining parameters.
 16. A method according to claim 15 wherein said laser machining device further comprises: an analysis unit of the state of said material to be machined connected to said central unit to provide a machining result, and that the following additional steps are implemented: e) acquiring with said analysis unit after step d), machining results; f) transmitting said machining results and said machining parameters to said central unit, said central unit being configured to determine information relating to the material to be machined; g) enriching said material database with said information relating to the material to be machined from said central unit. 17-18. (canceled)
 19. A laser machining system comprising means adapted to perform the steps of the method according to claim
 1. 20. The system according to claim 19 further comprising an optical unit for directing said laser beam towards said material to be machined.
 21. The system according to claim 20 wherein said optical unit for directing said laser beam towards said material to be machined allows the precession of said laser beam.
 22. A computer program for implementing a method according to claim
 1. 23. (canceled)
 24. A computer-readable medium on which the computer program according to claim 22 is recorded. 